Measuring bipartite quantum correlations of an unknown state 
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We report the experimental measurement of bipartite quantum correlations of an unknown two-qubit state. 
Using a liquid state Nuclear Magnetic Resonance (NMR) setup and employing geometric discord, we evaluate 
the quantum correlations of a state without resorting to prior knowledge of its density matrix. The method is 
applicable to any 20d system and provides, in terms of number of measurements required, an advantage over full 
state tomography scaling with the dimension d of the unmeasured subsystem. The negativity of quantumness is 
measured as well for reference. We also observe the phenomenon of sudden transition of quantum correlations 
when local phase and amplitude damping channels are applied to the state. 
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Introduction. — Quantum Mechanics rules the physical 
world, but sometimes its effects can be elusive. In particu- 
lar, entanglement, the genuinely quantum type of correlation 
shared by two parts of a composite system, is not directly mea- 
surable in laboratory. Indeed, there is not a self-adjoint opera- 
tor quantifying the amount of entanglement of a state. Hence, 
an a priori knowledge of the density matrix appears neces- 
sary to evaluate entanglement. This is a serious drawback, 
since dealing with high dimensional systems makes the state 
reconstruction extremely demanding if not unfeasible in terms 
of required resources. The problem has been overcome by in- 
troducing nontrivial lower bounds to entanglement measures, 
expressed as nonlinear functions of the density matrix coef- 
ficients, whose values can be detected by means of a limited 
number of measurements ||Tl|2l. 

However, it has been acknowledged that even separable 
(unentangled) mixed states show a distinctive nonclassical be- 
havior, benchmarked by their quantum discord |I3] |4l. Ar- 
guably considered the most general form of quantum correla- 
tions, it quantifies the minimum induced disturbance on a state 
of a bipartite system when a local measurement is made on 
one subsystem. Currently, massive theoretical and experimen- 
tal efforts are dedicated to discover the potential of such corre- 
lations for quantum information processing ll5 HT4l . The rea- 
son is that discord has appealing properties: For pure states, 
it reduces to entanglement; while for mixed states, even un- 
entangled states (apart from a subset of null measure) have 
nonvanishing discord fT5\. A number of phenomenological 
investigations have proven discord to be easier to create and 
far more robust than entanglement under decoherent dynamics 
116 1 . In spite of that, even if at theoretical level a plethora of 
discord-like measures have been introduced [5^ JJ|, and some 
discord witnesses have been implemented lfT8l420l . the situ- 
ation is similar to entanglement: There is not an observable 
quantity grasping the amount of such general quantum corre- 
lations in a state. One could try to recast a measure as a func- 
tion of density matrix elements and therefore of observable 



quantities, but most known quantifiers are entropic, thus can- 
not be easily associated to observables. Taking into account 
distance-based measures offers a solution to such conundrum. 

In this paper, following the theoretical proposal of 
Ref. Il2n . we report the experimental measurement of quan- 
tum correlations of unknown two-qubit states in a room- 
temperature NMR system Il22l - E5l . In this setting, the infor- 
mation is stored in the nuclear spins, while transformations 
and state preparation are implemented by applying highly 
controllable radiofrequency (rf) pulses, magnetic field gradi- 
ents and evolutions under spin interactions ll26l . On the other 
hand, the environment affects the spin system by inducing re- 
laxation that drives the system back to the thermal equilibrium 
distribution according to two independent characteristic times 
Ti, T2 El- 

We pick primarily the geometric discord Dq as quantifier 
of bipartite quantum correlations L28J . Its peculiarity is that 
it can be expressed, for arbitrary states of 2® d systems, as 
a state-independent function of the density matrix elements, 
and consequently in terms of observables |21, 29). In fact, 
there is no need to have full information on the state in order 
to measure geometric discord. We experimentally verify the 
theoretical results by preparing an unknown two-qubit Bell 
diagonal state (the choice of the state is due to experimen- 
tal convenience exclusively) and retrieving the value of geo- 
metric discord by means of local measurements over one of 
the subsystems. Then, we investigate the robustness of quan- 
tum correlations under phase damping and amplitude damp- 
ing channels acting separately on each qubit. It is predicted 
that, by appropriately engineering the initial state, measures of 
discord should undergo a sudden transition in their dynamical 
evolution 1 30 , 3 1 1 , exhibiting different regimes of resilience to 
decoherence under noisy conditions. We carry out a compre- 
hensive analysis of general quantum correlations. We monitor 
in fact the evolution of geometric discord ||28l , of its lower 
bound defined in ETIl . as well as of the so-called negativity of 
quantumness l32\ for the produced states (the latter measure. 
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here investigated for the first time in open systems), linking 
the singularity in their expression with the discontinuity in the 
dynamics. 

It is important to discuss the consistency of geometric dis- 
cord Dc as a reliable estimator of quantum correlations. Re- 
cent works have shown that measures built on the Hilbert- 
Schmidt norm, as geometric discord is, can increase under 
local reversible operations on the unmeasured subsystem, be- 
ing biased by the purity of the global state Il33l[34l . In spite of 
that, geometric discord is still a useful signature of quantum 
correlations in a number of relevant cases. For example, it can 
be legitimately used to investigate correlations between sys- 
tem and environment in open quantum evolutions |35 1, which 
are globally unitary thus leaving the purity unchanged. Even 
for a bipartite state under local decoherent evolutions, as stud- 
ied in this paper, the geometric discord reliably identifies the 
sudden transition point in the dynamics as well as other full 
fledged measures of discord ||3] [30l |3T| . Specifically, for a 
system of two qubits (keeping the dimension fixed), as in our 
case, the problem highlighted in |33| cannot occur, and Dc 
admits moreover an operational interpretation in terms of fi- 
delity of remote state preparation |10|, while also bounding 
from above 1,36,1 the entanglement quantified by the negativity 
||37l . In general, Dq is a valid lower bound to measures of 
quantum correlations based on relative entropy |5 1. 

As anticipated, we also investigate experimentally two 
other discord-like quantities. First, we take in account a non- 
trivial experimentally-friendly lower bound Q to geometric 
discord [21 J: such a quantity is less accurate and does not 
reveal any of the mentioned dynamical phenomena, as its ex- 
pression is continuous |38|. Then, we consider as a reference 
a fully bona fide measure of discord, the negativity of quan- 
tumness Q'^ introduced in 191 and discussed in detail in ll32l . 
Such a measure quantifies the minimum negativity |37l cre- 
ated with an apparatus during a local measurement over one 
subsystem of a bipartite system. It can be alternatively inter- 
preted (when the measured subsystem is a qubit) in terms of 
minimum trace distance from the set of zero-discord states. 
The Qj^ is unfortunately harder to compute and less acces- 
sible experimentally; to measure it from direct data, without 
tomography, we need a partial knowledge of the form of the 
state, specifically the fact that it is a Bell diagonal state. Such 
a measure is found to detect the same dynamical features as 
geometric and entropic discord (see also |39|). 

Theoretical predictions. — In the following, we consider a 
bipartite state of a 2 ® t/ system, described by a density matrix 
written in the Bloch representation [40]: 
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tions, while x - {xy},y - {y^} are the column vectors asso- 
ciated to each subsystem, and C = (Cy^) is the correlation 
matrix. Quantum discord evaluates the minimum induced 
disturbance on the state by a measurement over one subsys- 
tem, say A. One can adopt a geometric point of view and 
quantify quantum correlations as the minimum distance of 
the state to the set of classical-quantum states, i.e., the zero- 
discord states, which take the form x - 2; Pi\i){i\A ® Pbi, 
where YjiPi - 1 ™d d')) is an orthonormal vector set ||5]. 
In this context, the geometric discord introduced in Ref. Il28ll 



is defined as Dcip) - 2min^||p - xW^ where the squared 
Hilbert-Schmidt norm \\M\\\ - Tr[M"'"M] is picked as met- 
ric. An explicit formula of geometric discord for 2 ® c/ states 
is Dg(p) = idil' + ||C||2 - Ak^,,) = 2(Tr[5] - ^^ax), where 
^max is the largest eigenvalue of the matrix 5 = ^(xr^H-CC^). 
Expressing A;max explicitly, the geometric discord can be recast 
as a state-independent function of the density matrix elements 
lEI: 

Dc = ^Tr[5] - ^ V6Tr[52] - 2Tr[5]2 cos (^) , (1) 

where 9 - arccos {V2(2Tr[5]^ - 9Tr[5]Tr[52] -i- 

9Tr[5^])(3Tr[52] - Tr[5]2)"'}. Fixing = 0, one ob- 
tains a tight lower bound Q to the geometric discord, 
which is itself a faithful estimator of quantum correlations 
ETIl . The result allows to recast Dc and Q as functions 
of the expectation values of a set of observables [O,]: 
Dc = .foaiiOi)], Q = foWi)] EHl. The choice of the specific 
operators depends on the experimental setting. A proposal for 
an all-optical setup has been advanced in Refs. 11211 |29ll . In 
the present work, an NMR system is used. In this context, the 
column vectors x and y are proportional to the magnetization 
of each nuclear spin, x - 2{I ® I2) and y - 2(l2 ® /), where 
/ - {I^,Iy,I,} is the nuclear spin operator, which for spins 
- 2 is / - cr/2. The elements of the correlation matrix 

CyA - (o"v®cr^) - 4(/v®/i) andi^= 2(/®l2) are proportional 
to a combination of multi quantum coherences [23, 24|. 

A full state reconstruction demands 4-d^ - 1 spin measure- 
ments. However, to compute Dc and Q one can get rid of 
d^-l local measurements on one the subsystems, as the Bloch 
vector y does not account for the S matrix. Indeed, quan- 
tum correlations are evaluated from observables of the type 
(OyA) = Tr[o-y®TAp],v = I,..., 3; A = 0, . . . ,d^ - I, hut 
global measurements can be replaced by local ones ll25l : 



where RyA = Tr [p(cry ® cr^)], ctq = h,o-v (v = 1,2,3) are 
the Pauli matrices and {ta} their ^/-dimensional generaliza- 



Tr[(cr, ® crA)p] = Tr[(cri ® Ij)^,,] 

^yA = UyApUl, (2) 

where UyA = Ka^bR^^,^MaI for R<p,.,<i,,iOvA) = R^i^vA) ® 
R^i^OyA), being R^^^^^(OyA) a rotation by an angle OyA over the 
direction (fiy^), the indexes v,A - 1,2,3 refer to the rotations 
for measuring the C matrix elements and Ka^b represents the 
CNOT gate with A being the control qubit (see the Supple- 
mental Material [43] for details). 

The negativity of quantumness is defined for 2 ® li sys- 
tems as the minimum trace distance from the set of classical- 
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quantum states (see ll32l for other definitions and interpreta- 
tions), e^(p) = i min^j, llp-;rlli. where ||M||i = Tr[ Vm^m] is 
the trace norm. For the particular case of Bell diagonal two- 
qubit states, Q'^ is analytically computable |[32l l4r|. Denoting 
the ordered singular values of the Bloch correlation matrix C 
as |c,| > \cj\ > IqI, where /, j, k are permutations of 1, 2, 3, then 
the negativity of quantumness is given by half the intermedi- 
ate one, - \cj\l2. Therefore the negativity of quantumness 
can be also experimentally determined as described above. 

Experimental results. — In our implementation we set d - 
2. The state of a NMR two-qubit system in the high tem- 
perature approximation is given by p - \Ia + sAp, where 
s = ha>i/4kBT ~ 10"^ is the ratio between the magnetic and 
thermal energies, a>L is the Larmor frequency, kg is the Boltz- 
mann constant and T the room temperature [22J. All measure- 
ments and transformations affect only the deviation matrix Ap, 
which in fact contains the information about the system state. 
The unitary operations over Ap are implemented by radiofre- 
quency pulses and evolutions under spin interactions with an 
excellent control of the rotation angle and direction. NMR 
probes the transverse magnetization, which is proportional to 
the average values of the (/ J or (/,) operators, so only a few 
elements (single quantum coherences) are directly accessible. 
In this direction, to obtain a full characterization of Ap it is 
necessary to execute a set of independent measurements after 
applying specific rotations to the system, which characterize 
NMR quantum state tomography 1,24 J . It is worth remarking 
that, since in NMR experiments only the deviation matrix is 
detected, the calculations of the x vector and of the correlation 
matrix C are done in units of e. 

The experiments were performed on a liquid state carbon- 
13 enriched chloroform sample (CHCI3) at room tempera- 
ture, with the two qubits being encoded in the 'H and '^C 
spin-j nuclei |43|. Two initial Bell diagonal states pi_2 - 
l(h + H, cP''*^V; ® o"/) were prepared by mapping the correla- 
tion matrix into the deviation matrix, as described in 1 (311 l43l . 
Note that the methodology here employed is not restricted to 
this class of states, but extends to arbitrary states of 2 ® li sys- 
tems. To verify the correct preparation of the states, the pulse 
sequences proposed in Refs. [19. .31 1 were applied followed 
by a quantum state tomography procedure ll42]| . The result- 
ing deviation matrices are shown as block diagrams in Fig. [T| 
The interaction with the environment in a NMR system is de- 
scribed by phase and amplitude damping channels B4l . As 
can be seen in Refs. fT9','3Tl, the correlations presented in the 
initial state of Fig{T|a) should decay monotonically in time, 
while the state of FigjTJb) is expected to exhibit a nontrivial 
behavior of quantum correlations during the evolution. 

According to Eq. ([TJ and the accompanying discussion, 
quantum correlations in a two-qubit system are obtained from 
the expectation values of the correlation matrix elements Cy^ - 
{cTy ® cTji} - Tr[(cry ® crOp] and x - {if ® I2). These can be 
calculated for instance from the experimentally tomographed 
deviation matrix. However, using the rotations described by 
Eq. (j2]) with proper set of angles, the evaluation of Dq and 
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FIG. 1. (color online). Block diagrams for the deviation density 
matrix related to the Bell diagonal initial states. For pi in FigfTJa) we 
have \cf\ = |4"l = l4'*l = 0-2, while p2 in Figjljb) reads |cf T= 0.5, 
\cf\ = 0.06, and \cf\ = 0.24. 



Q is reduced to a set of spin magnetization measurements on 
one of the qubits. This means that, after a correct set of ro- 
tations (rf pulses) applied to the prepared state, the geometric 
discord (and its lower bound) can be determined directly from 
the NMR signals (apart from the factor), without having to 
know the state. We dub this procedure as direct measurement. 
The procedure is repeated until the thermal equilibrium state 
is re-established. Concerning the negativity of quantumness 
II32II . a partial knowledge of the state is instead required to 
evaluate its value. In particular, assuming (as a posteriori ver- 
ified by tomography, see |43|) that the state remains in Bell 
diagonal form during the evolution, then the direct method 
still suffices to extract the correct value of Q^. To benchmark 
the effectiveness of the direct method, overall, we compared 
its outcomes with the results obtained for the corresponding 
measures of quantum correlations by evaluating them on the 
state reconstructed by complete tomography as well. 

The results obtained for Dc and Q using the direct mea- 
surement and tomography procedures as well as the theoreti- 
cal predictions are reported in Fig. |2] for both experimentally 
produced states. In Fig|3] we present the results obtained for 
by direct measurements (the '// nucleus was detected), 
tomography and theoretical predictions for both states. We 
highlight a satisfactory agreement between the direct mea- 
surements, the tomographic data, and the theoretical predic- 
tions. This demonstrates that we are able to directly quantify 
bipartite quantum coiTelations in unknown (or partially known 
in the case of Q^) two-qubit states with our NMR setup. In 
particular, for the first time the negativity of quantumness is 
observed to undergo sudden transition in the same dynamical 
conditions as the geometric discord and the entropic discord 
Il30ll39l . Furthermore, directly assessing the value of Dc pro- 
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FIG. 2. (color online). Time evolution of Q and Dq for the state pi in 
(a) and (b), and p2 in (c) and (d). The green dots represent the direct 
measurement and the black open squares corresponds to tomography 
results. The red lines depict the theoretical predictions. The relax- 
ation times retrieved from experimental data are: Ti = 3.57 s and 
T2 = 1.2 ^ for the hydrogen, and Ti = 10 s and T2 = 0.19 i for the 
carbon. Quantum correlations are displayed in units of er. 



vides a nontomographic method to estimate a priori the suit- 
ability of the produced states for quantum communication via 
remote state preparation ifTOl . 

Conclusions. — NMR systems are a natural arena for quan- 
tum information processing with negligible entanglement, 
thus ideal testbeds for investigating dynamical properties of 
quantum correlations in open system dynamics. In this re- 
spect, the information stored in the density matiix of a state 
is often redundant to the quantification of its quantum corre- 
lations. Recent theoretical findings allow one to implement 
optimized protocols by reducing the number of required mea- 
surements as compared to the full state reconstruction. Here 
we have considered an NMR system and detected the quan- 
tum correlations of a two-qubit Bell diagonal state. Further- 
more, we observed the sudden transition of geometric discord 
and negativity of quantumness under phase and amplitude 
damping channels. The sudden transition and in particular the 
freezing ll30ll . common to various measures of quantum cor- 
relations other than entanglement under particular decoherent 
evolutions, is certainly worthy of further theoretical and ex- 
perimental investigation 1451 . 

We finally stress that the reported experiment has to be 
taken as a proof of principle. Thanks to the explicit closed for- 
mula, without any knowledge of the initial state of a 2®(i sys- 
tem, it is possible to evaluate geometric discord by 3d^ local 
spin measurements, and tune the dynamics in order to protect 
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FIG. 3. (color online). Time evolution of |ci| (upward triangle), |c2| 
(downward triangle) and IC3I (hexagon) for state pi in (a), and p2 
in (c). The open symbols represent the tomography results and the 
filled symbols represent the direct measurement. Note that for each 
|c,| (i = 1, 2, 3) both forms of measurement are superimposed. In (b), 
for the initial state pi, and (d), for the initial state p2, the green dots 
represent the direct measurement and the black squares represent the 
tomography results. The red lines are the theoretical predictions for 
Q^. Quantum correlations are displayed in units of s. 



quantum correlations from decoherence. Also, it is remark- 
able that an optical setup would allow to detect quantumness 
of correlations by only seven projective measurements for any 
d [2T\, but the implementation of such quantum operations ap- 
pears challenging even with state-of-the-art technology |38|. 
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SUPPLEMENTAL MATERIAL 
NMR experiments 

The NMR experiments were performed on a carbon- 13 en- 
riched chloroform sample (CHCI3), where the two nuclear 
spins of 'H and '^C implement the qubits. The sample was 
prepared by mixing 100 mg of 99 % '^C-labelled CHCI3 in 0.7 
mL of 99.8% CDCI3 in a 5 mm NMR tube (both compounds 
provided by the Cambridge Isotope Laboratories Inc.). NMR 
experiments were carried out at 25°C in a Varian 500 MHz 
Premium Shielded spectrometer located at the Brazilian Cen- 
ter for Research in Physics (CBPF, Rio de Janeiro). A Varian 
5 mm double resonance probehead equipped with a magnetic 
field gradient coil was used. Two njl pulses of length 7.4 /is 
and 8.0 jUS were applied to 'H and '^C, respectively. Spin- 
lattice relaxation times (Ti) for 'H and '^C nuclei, measured 
by the inversion recovery pulse sequence, were 3.57 s and 
10 s. Spin-spin relaxation times (T2), measured by a CPMG 
pulse sequence, were estimated to be about 1.2 s and 0.19 s 
for 'H and '^C nuclei. The recycle delay (di) was set at 60 s 
in all experiments. 

The nuclear spin Hamiltonian is given by 

'H^-iojH- w^) /f - (wc - wf^) + 2nJ /f /f (A.3) 
-I- {l" cos ip" + l" sin ip") + (/f cos p'^ + if sin cp^) , 

where ^/^) is the spin angular momentum operator in the 
a,j3 = jc, y, z direction for 'H ('^C); ip'^ (ip''^ defines the direc- 
tion of the rf field (pulse phase) and ('^f) the rf nutation 
frequency (RF power) for 'H ('^C) nuclei. The first two terms 
describe the Zeeman interaction between the 'H and '^C nu- 
clear spins and the main magnetic field Bq. The correspond- 
ing frequencies are conl^n ^ 500 MHz and <jJcl2n k 125 
MHz. The third term is due to a scalar spin-spin coupling 
of y » 215.1 Hz. The fourth and fifth terms represent the rf 
field to be applied to the 'H and '^C nuclear spins, respec- 
tively. There is also a time-dependent coupling of the nuclear 
spins with the environment that includes all fluctuating NMR 
interactions accounting for the spin relaxation, e.g., 'H - '^C 
dipolar spin-spin couplings and interactions with the chlorine 
nuclei. 

To investigate the dynamics of quantum correlations, we 
prepared a Bell-diagonal state: 



I +Ci ci - C2 ^ 

1 - C3 Ci -H C2 

Ci -H C2 1 - C3 

\Cx - C2 1 -H C3 



(A.4) 



where the coefficients {c,) are picked as explicited in Fig. 2(a) 
of the main text. Such state is obtained from thermal equi- 
librium by applying the pulse sequence for producing the 
pseudo-pure state |11), followed by rf pulses which imple- 
ment a pseudo EPR gate liAlil . It is worth mentioning that 



a true EPR gate would require additional pulses. To pre- 
pare the state showed in Fig. 2(b), a diagonal state was created 
from the thermal equilibrium by 5.5 ms strongly modulated 
pulses (SMP) - designed using a MATLAB self-written rou- 
tine based on a SIMPLEX optimization protocol - followed 
by a 2 ms magnetic field gradient pulse. The deviation ma- 
trix is Ap = diag(Q',y6, 7, 6), whose entries are defined through 
the optimized SMP pulse. Next, the deviation matrix is trans- 
formed by a pseudo-EPR gate into an X-type matrix MAll : 



a + y —a + y ' 

(3 + 6-/3 + 6 

-13 + 6 13 + 6 

—a + y a + y 



(A.5) 



The SMP was designed to produce a diagonal state with pop- 
ulations a, /3, y and 6, so that, after the transformation by the 
pseudo-EPR gate, a state with a deviation matrix in the form 
equivalent to the state in Eq. (2) of the main text, with spec- 
ified c,'s, was obtained. After having prepared the desired 
state, the system was left to evolve for a period t - m/4J, with 
m = 1,2,3,4,... The dynamics, considering on-resonance 
evolution, was characterized by two terms of the Hamiltonian 
shown in Eq. (A.4i acting on the nuclear spins: the scalar 



coupling and the time-dependent term, which account for the 
system relaxation. 

After the evolution period, to read-out the quantum state, 
we implemented the protocol of Fig|4] The y rotation was 
carried out by a single rf pulse, while for the z rotation the 
pulse sequence (|) was applied to both nuclei. 

The controlled-NOT gate was achieved by the pulse sequence 

(f):-"(iMft-(!£,-(!t-(f):-(f)",,-(f):-(f):.»'>- 

the super indices states for nucleus where the pulse is applied 
and U (^jj^ represent a free evolution under J coupling. In this 
case, a pulse sequence implementing a true CNOT gate is in- 
deed necessary to produce the correct output. As described 
above, the experimental procedure depicted in Fig. 1 was run 
three times for each initial state in order to measure the mag- 
netization (/f in the states ^, that leads to determine the 



two-point correlation functions {ij^ 



'If), 



So, the quantum 



correlations measured by Dq, Q, could be directly evalu- 
ated. The rotations to transform all of the elements of the cor- 
relation matrix into magnetization ones {{if ® lf}p — » (/f )f,) 
are given by: 





= RAO) 




*,(fe) 


= Rzzin/2) 




63(^33) 






62(^12) 


= R„(37t/2) 


R<I>1,<, 




= ^,,(37r/2) 


^01 ,( 




= R.y(7r/2) 


^1*3 ,( 


*,(03l) 


= Ryj7r/2) 




= -R,y(n/2) 


^03, ( 




= -Ry,i7T/2) 
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FIG. 4. Schematic representation of the pulse sequence for the direct 
magnetization measurements. The unfilled red bars represent pulses 
whose phases and angles depend on the indexes v, A = 1, 2, 3, which 
refer to the measure of each element of the C matrix. The phases 
and angles are 0i = x, 02 = z, 'p3 = y, Sii =0, 622 = 833 = ^/2, 
On - ^21 = 3;r/2, 613 = 631 = 7r/2 = 6*23 = ^32 = ^/2. The thinner 
blue bars represent a ;r/2 pulse realizing the CNOT logical gate. 



The quantum state was also read-out by quantum state to- 
mography, as described in Ref. MA2IIA31 . in order to check the 
consistency of the results obtained by direct measurements. 
This process allowed us to reconstruct the deviation density 
matrix of the system for each instant of time, benchmarking its 
behavior under the coupling with the environment. Multi-step 
increments of t in successive experiments allowed the effect 
of the spin environment on the initial deviation matrix to be 
followed. Final states corresponding to 25 1 distinct t values 
were acquired. The high homogeneity of the static magnetic 
field was guaranteed by the linewidth of * 0.95 Hz in 'H spec- 
tra. 



NMR presents, typically, two relaxation channels, namely 
the generalized amplitude damping and the phase damping. 

The amplitude damping channel is described by the follow- 
ing set of operators 

where, in the NMR context, y = 1/2 - e/2 and p = I - 
exp(-t/Ti), Ti being the longitudinal relaxation time of the 
qubit under consideration. We observe that in our case, the 
relaxation times are different for the two qubits, since they 
have distinct Larmor frequencies. 

For the phase damping channel, one has 

where A - 1 - exp {-tlTi) with Ti being the transverse relax- 
ation time associated with the qubit. 

The agreement between the model and the experimental 
data is quite satisfactory. In the theoretical model we ne- 
glected a number of possible error sources, such as mag- 
netic field inhomogeneities or oscillations of the residual off- 
diagonal elements of the deviation matrix. In spite of that, all 
the singular features predicted by the quantum model emerge 
from the experimental data analysis. Specifically, for Bell- 
diagonal states one has S - ^ - d\?ig{c\lA,c^lA,c^j^lA] and 
expressions for Dq, Q, are easily obtained. While the 
quantity Q is manifestly continuous, Dc and show dis- 
continuity in the quantum coiTelations at a critical time. 



Decoherence analysis 

The decoherence process is theoretically analysed through 
the operator sum representation technique |A4|, in which the 
evolution of the density operator is given by IA5I 

£>{PAB) = _^(£f ® ® £f )Pab(Ia ® £f )"'"(£f ® Ib)\ 
'J 

where E"^ (f) are the Kraus operators. 
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